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Abstract 



On the basis of a quantum-statistical approach to the electrical conductivity of nonideal 
plasmas we derive analytical results in the classical low-density regime, in the degenerate 
Born limit, and for the contribution of the Debye-Onsager relaxation effect. These explicit 
results are used to construct an improved interpolation formula of the electrical conductivity 
valid in a wide range of temperature and density which allows to compare with available 
experimental data of nonideal plasmas. 



Optical and transport properties of plasmas are governed by the mutual Coulomb 
interaction and depend strongly on temperature and the electron and ion density. Ex- 
perimental efforts for the measurement of the electrical conductivity up to high den- 
sities revealed the importance of collective excitations and many-body effects such as 
Pauli blocking, dynamic screening and self-energy, structure factor, the Debye-Onsager 
relaxation effect, formation of bound states etc. Although different methods have 
been proposed for the evaluation of plasma transport properties [fj, [| f|, ||, |J, it still 
remains a theoretical challenge to treat these effects by a unified quantum-statistical 
approach. 

On the other hand, in addition to highly sophisticated approaches to the conductiv- 
ity accounting for quantum statistical many-particle effects, for practical use interpola- 
tion formulas are of interest which are applicable in a large region of plasma parameter 
values. Such interpolation formulas were developed by several authors 0, ||, (|, §] in 
order to evaluate complex physical situations, e.g., in hydrodynamical simulations for 
the generation and expansion of plasmas produced by high-power lasers, energetic 
heavy ion beams or magnetic compression (pinches). Based on rigorous results for 
special limiting cases and possibly numerical simulations in the intermediate regions, 
Pade-type interpolation formulas have been proven to be highly effective to describe 
thermodynamic as well as tranport properties. 

In this paper we consider a fully ionized plasma consisting of electrons (mass m e ) 
and singly charged ions (mass M), interacting via the Coulomb potential, at density 
n e = rii = n and temperature T (hydrogen plasma). The dimensionless parameters 



1 Introduction 
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describe the ratio between the mean potential energy and the kinetic energy (r) and 
denote the degree of degeneracy of the electron gas (O). Usually, plasmas are clas- 
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sified into ideal (T <C 1) and nonideal (T > 1) as well as degenerate (0 -C 1) and 
nondegenerate (9 ^> 1) systems. 

Using these dimensionless plasma parameters, the construction of interpolation 
formulas for the dc conductivity has been performed in different works. Ichimaru 
and Tanaka |2| considered a two-component plasma at finite temperatures within a 
generalized Ziman formula by taking into account dynamic local field effects in the 
dielectric function and the dynamic ion-ion structure factor. The collision integrals 
were evaluated in Born approximation to study a strongly coupled and degenerate 
plasma with T > 1 and O < 1. The correct nondegenerate limit was incorporated 
in a semiempirical way by adopting a prefactor in correspondence to the well-known 
Spitzer formula. 

Ebeling et al. || and Lee and More [Q used the relaxation time approximation for 
the evaluation of the collision integrals both in the low-density (nondegenerate) and 
high-density (degenerate) regime. Electron-electron interactions are neglected in that 
approximation which is only justified in the degenerate case because of Pauli blocking. 
The consideration of electron-electron scattering is, however, decisive to obtain the 
right low-density result, the prefactor of the Spitzer formula. This has been modeled 
using a semiempirical ansatz. Furthermore, Rinker 1| has derived an interpolation 
formula for the electrical conductivity from an improved Ziman formula which are 
applicable in the strongly coupled, degenerate domain. 

Within a quantum statistical approach, the transport properties within the gener- 
alized linear response theory of Zubarev have been evaluated [ff| and rigorous analytical 
results for the electrical conductivity have been derived which are valid in the low- 
density (nondegenerate) and high-density (degenerate) limit J7|, ||, [|. The influence of 
non-equilibrium two-particle correlations in lowest order (Debye-Onsager effect) has 
been reconsidered recently [|ItJ. Based on these results and taking into account other 
approaches, we will construct an improved interpolation formula for the electrical con- 
ductivity to cover a wide range of temperatures and densities. 

2 Electrical conductivity 

Using the plasma parameters (|l|), we can express the electrical conductivity as 

(fc B r) 3 / 2 (4 7re 
m e e 



* {n ,T)= v^± T £p> „•(!■, 9) (2) 



with a universal function <r*(r, 0) depending on the characteristic plasma parameters 
exclusively. 

Our quantum-statistical approach to the electrical conductivity is based on a gen- 
eralized Boltzmann equation which is derived from linear response theory M. A finite- 
moment approximation of single and two-particle distribution functions defines the set 
of relevant observables B n , and the electrical conductivity is obtained from the ratio 



(3 1 

a = 



£2 \D[B n i] B n ] 



N[B n ] 
Q[B n ] D[B n ,-B n ) 



(3) 



The determinants contain equilibrium correlation functions defined by 

N[B n ] = (R; B n ) , (4) 
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Q[B n ] = (R; B n ) + (R; B n ) , 
D[B n i\ B n ] = {B n i\ B n ) + (B n >\ B n ) , 

which are given by Kubo's scalar product between two operators: 

(A;B) = \ [ drTr[p A(-ihr)B), (5) 
P Jo 



(A{e);B) = lim / dt e et (A(t); B) . 

e-*0 J-oo 



The two-particle center of mass momentum is defined via P = —em e M/(m e + M)R. 
Q is the system volume and (3 = 1/ksT. Furthermore, the time- dependence of the 
relevant observables B n and of the equilibrium statistical operator p is given by the 
standard system Hamiltonian of a fully ionized plasma, see Ref. J7j. The correlation 
functions D[B n r] B n ] are related to thermodynamic Green functions. Thus, efficient di- 
agram techniques provide a systematic treatment of many-particle effects in a strongly 
coupled plasma which has been demonstrated for the single and two-particle level, see 
Refs. [0,0. 

As the relevant observables on the single-particle level are treated the moments of 
the electron distribution (adiabatic approximation) 

P m = J2 hk ^{^) a H k ) a e(k), (6) 

including the total electron momentum (P ) and the ideal energy current (-Pi). A 
systematic increase in the number of moments results in a converging expression of 
the conductivity as it is known from the Kohler variational principle or the Grad and 
Chapman-Enskog method of kinetic theory. 

The correlation functions D[P n , P m ] are given by a sum of ladder diagrams in the 
low-density limit. Both the electron-electron and the electron- ion scattering processes 
are included, see Ref. ||. The scattering integrals have been evaluated in various 
approximations which are appropriate for a given density and temperature range. 

(i) The correlation functions in the classical limit (T 2 ^> 1) are related to Boltz- 
mann collision integrals which are given by transport cross sections. These quantities 
are treated in the quasi-classical approximation as relevant for low-energy particles. 
Additional quantum corrections are found from a WKB expansion of the collison inte- 
grals. Using a three-moment approximation for the one-particle distribution function 
[0 we find for T 2 9 > 1 

a*(r 2 0»l) = 2a (lnr- 3 + 0.2943 (7) 

The prefactor a = 0.591 is the limiting value of the moment expansion and coincides 
with the Spitzer result [fL^l - 

(ii) Born limit (T 2 1): The correlation functions are equivalent to the Lenard- 
Balescu collision integral if the random phase approximation (RPA) is considered 
for the dielectric function e^pA^q, lj). The following low-density limit is obtained for 
r 2 e < 1 0: 



a*(r 2 < 1) = 2a (in ® + 0.4807^) . (8) 
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(iii) The correlation functions in the high-density (liquid-like) limit are given by 
Landau collision integrals which treats the electron-ion scattering in Born approxima- 
tion. Thus, the corresponding electron-ion contribution is 



2u2\ m + n r2k 



Vei(q) 



2 



.(9) 



The number of required moments P m reduces with increasing density; the Ziman-Faber 
formula is already obtained from Eq. (^) with the lowest moment P . It includes the 
static ionic structure factor Su(q), the electron- ion pseudopotential V e i(q), and the 
static dielectric function of the electrons e(q). 

The consideration of two-particle correlations described by respective moments 

Sr C m \i) = / c m d m '( k > p> q)40 - \) a \i.p + §K(p - §K(fc + ~) (io) 

k,p 

leads to a decreased plasma conductivity. Spatial symmetry properties of the two- 
particle distribution function in the case of an applied electric field restricts the mo- 
ments to the electron- ion contribution 5n™' m (q) [|T0| | . The evaluation of the first 



moment f®f = 1 in the low-density regime gives the original Onsager result of the 
Debye-Onsager relaxation effect ]7[ pT0| . This effect decreases the conductivity due to 
an incomplete formation of the screening cloud and was first considered in the theory 



of electrolytes [14]. An equivalent formulation is the hydrodynamic approximation 



in kinetic theory and assumes a local equilibrium for the distribution functions [|T^ 



As this assumption already fails to describe the single-particle properties (Spitzer for- 
mula) properly also higher two-particle moments have been considered in Ref. [|l(| . 



Corrections to the Onsager result were indeed found for a fully ionized plasma from 
the second moments Sn1f(q) with f^f = k 2 and a virial expansion of the conductivity 

a* ^(lnr-^ + fr + cr^inr- 3 / 2 )- 1 (11) 

could be derived in the classical low-density and nondegenerate limit. We find for 
hydrogen plasmas 



0.591 



(47re ) 2 (fc B r) 

e 2 m l/2 



3/2 



lnr^ 3/2 + 1.124 + 



1.078 



p3/2 l n p-3/2 



(12) 



= 1.530 x 1(T 2 T 3/2 [lnr~ 3/2 + 1.124 + 0.258r 3/2 lnr~ 3/2 ] 1 (OmK 3 / 2 )" 1 . 

The inclusion of the two-particle nonequilibrium correlations determines the coefficient 
c of the term r 3 / 2 lnr" 3/2 in Eq. (0). 



3 Interpolation formula and results 

An interpolation formula can be constructed on the basis of the explicit analytical 
results given above so that reliable results for the electrical conductivity are obtained 
easily for a wide range of density and temperature without the necessity to evaluate 
the underlying complicated many-particle methods in full detail. The original inter- 
polation formula given in Ref. [S] included the low-density (f < 1) and Born limit 
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(r 2 B <C 1). We improve the validity region and the accuracy of that formula in the 
present paper by (i) taking into account the corrections due to the Debye-Onsager 
relaxation effect and (ii) by a better analytical structure of the interpolation formula, 
avoiding unphysical behavior such as, e.g., poles in the entire (T, 0)-plane. (iii) Fur- 
thermore we incorporate the results of Ichimaru and Tanaka in the strongly coupled 
and degenerate limit (0 < 1,T > 1) by analyzing their parameterized numerical re- 
sults. Thus the validity range of our interpolation formula is extended to a parameter 
range where the influence of the ion-ion structure factor becomes relevant. 
We propose the following interpolation formula [T in K, a in (Om) -1 ]: 



a T 3 / 2 ( 1 ' hl 



° = Q0T/ l 1 + 93/2 



m(l + A + B)D — C 



& 2 + re 



-i 



(13) 



with the functions 

A = r~ 3 

B 

C ■ 

D : 



l + a 2 /T^ + a 3 /T^ K + + 1)] , (14) 

ft 3 (i + c 3 e) 
re(i + c 3 e 4 / 5 ) ' 

in(i + r- 1 ) + c 5 r 2 e ' 
r 3 + a 5 (i + a 6 r 3 / 2 ) 
r 3 + a 5 



The parameters a, are fixed by the low-density virial expansion. In particular, the 
corrections from the Debye-Onsager relaxation effect are included in the function D. 
The coefficients 6j are used to adjust the Ichimaru and Tanaka results in the strong 
degenerate limit, and the parameters c$ were fitted to numerical data of the correlation 
functions. The explicit set of parameters is given by ao = 0.03064, ai = 1.1590, 
a 2 = 0.698, a 3 = 0.4876, a 4 = 0.1748, a 5 = 0.1, a 6 = 0.258, b x = 1.95, b 2 = 2.88, 
63 = 3.6, ci = 1.5, C2 = 6.2, C3 = 0.3, C4 = 0.35, C5 = 0.1. 

We compare the results of the new interpolation formula (0) with the former one 
given in 0, the Ichimaru- Tanaka fit formula 0, and experimental data for strongly 
coupled plasmas |1^, [T?|, |18| in Table [l|. Taking into account an experimental error 



of about 30% we find a good agreement between theory and experiment. However, 
complete correspondence with the experiments can not be anticipated because of de- 
viations from the Coulomb potential for the electron-ion interaction in rare gases and 
the occurrence of neutral particles (partial ionization) not included in the present 
theoretical approach, which is focussed to the fully ionized Coulomb system. The 
Debye-Onsager relaxation effect leads to a reduction of the electrical conductivity in 
the order of 5% at T w 1. A direct comparison with simulation results would be 
highly desirable. However, the standard approach through the current-current auto- 
correlation function is still limited by the small number of simulated particles and the 
corresponding statistical uncertainties. 

The validity range of the new interpolation formula (0) is restricted to a param- 
eter range where the formation of bound states (O > 1) can be neglected. Partial 
ionization plays a crucial role in low-temperature plasmas and can lead to a minimum 



in the isotherms of the electrical conductivity, see Ref. [19| for the case of hydrogen 



Bound states (atoms) are also important when evaluating the hopping conductivity 



6 



Contrib. Plasma Phys. () 



Tab. 1: Electrical conductivity according to the new interpolation formula ( J13| ) compared 
with the former one given in P], the Ichimaru-Tanaka fit formula 0], and available experi- 
mental data for strongly coupled plasmas 17, 18|. 



n, 



e 



Aril 



Xe 16 



Ne 16 



Air [H] 



H 18 



;i0 3 K) (10 25 /m 



a 

exp. 



TcFTr 

[2] 



[9] (13) 



22.2 
20.3 
19.3 
19.0 
17.8 
30.1 
27.0 
25.1 
22.7 
19.8 
19.6 
11.0 



Ar 17 16.4 



Xe|7| 12.4 
12.6 



15.4 
18.7 
21.5 



2.8 

5.5 

8.1 

14 

17 

25 

79 

160 

200 

1.1 

1.9 

0.13 



0.06 

0.1 

0.13 

0.15 

0.06 

0.12 

0.07 

0.14 



0.1 

0.15 

0.25 



0.368 

0.505 

0.604 

0.736 

0.838 

0.564 

0.922 

1.26 

1.50 

0.303 

0.367 

0.267 



56.9 
33.2 
24.4 
16.7 
13.7 
17.9 
7.47 
4.34 
3.38 
94.6 
65.0 
218 



190 
155 
170 
255 
245 
450 
740 
780 
930 
130 
165 
60 



200 
203 
209 
234 
232 
442 
546 
660 
694 
148 
160 
53 



225 
231 
239 
269 
270 
458 
594 
813 
933 
173 
187 
67 



0.128 

0.165 

0.18 

0.19 

0.185 

0.234 

0.192 

0.239 



551 
385 
324 
291 
403 
252 
371 
238 



83 
79 
76 
64 
46 
41 
48 
44 



78 
83 
86 
88 
55 
60 
57 
65 



93 
102 
105 
108 
70 
76 
73 
79 



0.175 
0.165 
0.170 



364 
337 
276 



62 
91 
114 



77 
103 
131 



95 
125 
156 



220 
224 
232 
261 
262 
453 
590 
797 
900 
169 
181 
65 



93 
102 
105 
107 
69 
75 
72 
77 



94 
124 
155 



for fluid hydrogen near to the nonmetal-to-metal transition at megabar pressures [^(J 



At present, the treatment of the dc conductivity for plasmas where bound states may 
occur is performed within the model of the partially ionized plasma, where the com- 
position is determined within a thermodynamic approach, and the contribution of 
the interactions between the different components is considered separately. Within 
such an approach, the interpolation formula given here may be of use to describe the 
contribution of the interaction between electrons and ions to the conductivity. 

The generalization of the present approach to plasmas with higher charged ions 
Z > 1 is also possible but not intended here. Then, comparison with new experimental 
data for strongly coupled metal plasmas with r 1 0, |22| can be performed, see 



23j. The inclusion of bound state formation and the extension to higher charged ions 
should be considered as possible goals of future work on interpolation formula for the 
electrical conductivity of nonideal plasmas. 
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